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Abstract 



RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to 
both monovalent and multivalent ions. Despite recent advances, we do not have general methods 
to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, 
and the thermodynamic parameters for secondary structure prediction have only been derived 
at 1M [Na + ]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical 
tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed 
monovalent /Mg 2+ salt conditions. We measure the unfolding and folding rupture forces and 
apply Kramers theory to extract accurate information about the hairpin free energy landscape 
under tension at a wide range of ionic conditions. We obtain non-specific corrections for the 
free energy of formation of the RNA hairpin and measure how the distance of the transition 
state to the folded state changes with force and ionic strength. We experimentally validate the 
Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test 
the approximate rule by which the non-specific binding affinity of divalent cations at a given 
concentration is equivalent to that of monovalent cations taken at 100 fold that concentration 
for small molecular constructs. 



1 Introduction 



RNA hairpins are elementary structures found in many macromolecular assemblies. It is generally 
accepted that a deeper understanding of their dynamics is a critical step towards the elucidation 
of many biological processes, like the regulation of gene expression [TJ [2 [2 HI E] ; the catalytic 
activity in many reactions jSJ [7] ; the ligand-binding specificity [5] ; or the RNA folding problem 
[HI HO] • DNA and RNA hairpins are also appealing model systems for their simplicity as they 
are amenable to exhaustive studies using a more physically-oriented approach, where theoretical 
models can be rigorously tested using simulations and experiments [111 112] . 

Many different and complementary biophysical methods have been used to study these struc- 
tures. For example, using time-resolved nuclear magnetic resonance (NMR) spectroscopy and 
thermal denaturation experiments, kinetics and thermodynamics of bistable RNA molecules were 
studied [131 [111 [111 [T71 HU HOI HH [H] . Recently, a photolabile caged RNA was designed 
to stabilize one ground-state conformation and study the folding kinetics by NMR and CD spec- 
troscopy under different conditions, including Mg 2+ |23j . Laser temperature-jump experiments 
have also been used to characterize the folding kinetics of small RNA hairpins at the ns and 
/is timescales [531 UHl HZ1 I2H1 US]- Using coarse-grained Go- like models, it was predicted 
that hairpins unfold in an all-or-none process in mechanical experiments [30], in agreement with 
experimental results [311 E21 [33] . 

Within the cell, many dynamical processes involving transient melting events of DNA and 
RNA double strands are driven by the application of localized forces by molecular motors. There- 
fore, single-molecule experiments are ideal to understand the thermodynamics and kinetics of 
macromolecules inside cells [3H [35] . As pointed out by Hyeon et al. .30J , force-denaturation us- 
ing single-molecule experiments are intrinsically different from thermally-induced denaturation: 
in bulk experiments where the unfolded state is accessed by raising the temperature or lowering 
the concentration of ions, the unfolded state is a high-entropy state while in mechanical pulling 
experiments the unfolding process is a transition from a low-entropy state to another low-entropy 
state. Regions of the free energy landscape normally inaccessible by conventional methods are 
probed using mechanical experiments. Consequently, pathways and rates of thermally-induced 
and mechanical unfolding processes are expected to be different. 

In a previous work [36] we pulled an RNA hairpin using optical tweezers [37l [38] to study the 
base-pairing thermodynamics, kinetics and mechanical properties at a fixed monovalent condi- 
tion. A kinetic analysis was introduced to determine the location of the force-dependent kinetic 
barrier, the attempt rate, and the free energy of formation of the molecule. Here we performed a 
systematic study by mechanically pulling the same RNA hairpin at different monovalent cation 
concentrations and also at mixed ionic conditions containing different concentrations of Mg 2+ 
cations. This is important because RNAs are highly charged polyanions whose stability strongly 
depends on solvent ionic conditions [39 [ 140" ! |4T [ |42" 1 143] . Despite its biological significance, we have 
limited information about RNA helix stability in mixed monovalent/multivalent ionic conditions 
[331 125] ■ In fact, the thermodynamic parameters for secondary structural elements of RNAs have 
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only been derived at the fixed standard salt condition of 1M [Na+] [TBI 115] IT? ] |Tg| \T§ \ ISO ] |2"T] . 
Here we derived numbers such as the persistence length describing the elastic response of ssRNA 
and also the free energy of formation of an RNA hairpin at different monovalent and mixed 
monovalent/Mg 2+ conditions. Our results are compatible with predictions obtained using the 
Tightly Bound Ion (TBI) model for mixed ion solutions, that treats monovalent ions as ionic 
background and multivalent ions as responsible from ion-ion correlation effects, and which takes 
into account only non-sequence-specific electrostatic effects of ions on RNA [15 } 144 } 135 ] . Our find- 
ings demonstrate the validity of the approximate rule by which the non-specific binding affinity 
of divalent cations is equal to that of monovalent cations taken around 100 fold concentration 
for small molecular constructs [JSJ H7] . 

2 Materials and Methods 

2.1 Molecular Synthesis 

The RNA molecule was prepared as previously described |3S]. Oligonucleotides CD4F (5'- 
AATTCACACG CGAGCCATAA TCTCATCTGG AAACAGATGAG ATTA TGGCTCGC ACACA- 
3') and CD4R (5'-AGCTTGTGT GCGAGCCATA ATCTCATC TGTTTCCAGAT GAGAT- 
TATGGC TCGCGTGTG-3') were annealed and cloned into the pBR322 DNA plasmid (Gen- 
Bank J01749) digested with EcoRI (position 4360) and Hindlll (position 30). The annealed 
oligonucleotides contain the sequence that codes for a modified version of CD4-42F class I hair- 
pin that targets the mRNA of the CD4 receptor of the human immunodeficiency virus [49] , 
Oligonucleotides T7_Forward (5'-TAATACGACTCA CTATAGG GACTGGTGA GTACTCA 
ACCAAGTC-3') and T7_Reverse (5'-TA GGAAGC AGCCCAGT AGTAGG-3') were used as 
primers to amplify by PCR a product of 1201 bp from the recombinant clone containing the 
CD4 insert. This amplicon contains the T7 RNA Polymerase promoter at one end, and was 
used as a template to synthesize an RNA containing the RNA hairpin (20bp stem sequence and 
tetraloop GAAA) and the RNA components of handles A (527 bp) and B (599 bp). The DNA 
components of handles A and B were obtained by PCR from the pBR322 vector (positions 3836-1 
for handle A and positions 31-629 for handle B). Handle A was 3' biotinylated while handle B 
was tagged with a 5' digoxigenin. Hybridization reactions were performed in a formamide-based 
buffer |5D] with a step-cool temperature program: denaturation at 85°C for 10 min, followed by 
1.5 h incubation at 62°C, 1.5 h incubation at 52 °C, and finished with a cooling to 10°C within 
10 min. 

2.2 Measurement Protocol 

All experiments were performed using a dual-beam force measuring optical trap [37l [38] at 25±1 
°C in buffers containing 100 mM Tris.HCl (pH 8.1), ImM EDTA, and NaCl concentrations 
of 0, 100, 500, and 1000 mM, or in buffers containing 100 mM Tris.HCl (pH 8.1) and MgCl 2 
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concentrations of 0.01, 0.1, 0.5, 1, 4, and 10 mM. The monovalent cation concentration [Mon + ] 
includes the contributions from [Na + ] ions and dissociated [Tris + ] ions. At 25°C and pH 8.1, 
about half of the Tris molecules are protonated, therefore 100 mM Tris buffer adds 50 mM to the 
total monovalent ion concentration |51j . Anti-digoxigenin polyclonal antibody-coated polystyrene 
microspheres (AD beads) of 3.0-3.4 /j,m (Spherotech, Libertyville, IL) were incubated at room 
temperature with the molecular construct for 20 min. The second attachment was achieved 
inside the microfluidics chamber using a single optically trapped AD bead previously incubated 
with the RNA hairpin and a streptavidin-coated polystyrene microsphere (SA bead) of 2.0-2.9 
/im (G. Kisker GbR, Products for Biotechnologie) positioned at the tip of a micropipette by 
suction (Fig. 1A and IB). 

Tethered molecules were repeatedly pulled at two constant loading rates of 1.8 pN/s or 12.5 
pN/s by moving up and down the optical trap along the vertical axis between fixed force limits 
and the resulting force-distance curves (FDCs) were recorded (Fig. 2A). A pulling cycle consists 
of an unfolding process and a folding process. In the unfolding process, the tethered molecule 
is stretched from the minimum value of force, typically in the range of 5-10 pN, where it is 
always at its native folded state, up to the maximum value of force, typically in the range of 
25-30 pN, where the molecule is always unfolded. In the folding process the molecule is released 
from the higher force limit (unfolded state) up to the lower force limit (native folded state) [55] . 
A minimum of two molecules (different bead pairs) were tested at each ionic condition, and a 
minimum of 100 cycles were recorded in each case (detailed statistics are given in the Supporting 
Material, section SI) 

2.3 Hairpin Model 

Under applied force it is feasible to reduce the configurational space of an RNA hairpin containing 
N base pairs (bps) to a minimum set of N + 1 partially unzipped RNA structures [3BJ [531 |53] . 
Each configuration in this set contains n adjacent opened bps in the beginning of the fork followed 
by N — n closed bps, with < n < N. The folded state (F) is defined as the configuration in 
which n = (all bps are formed), and the unfolded state (U) is the hairpin configuration in which 
n = N (all bps are dissociated) . Based on a simple calculation (see Supporting Material, section 
S2) we conclude that fraying |55j plays a rather minor role (if any) on the folding/unfolding 
kinetics of the sequence under study (Fig. 1A) and we do not include it in our analysis. The 
stability of each configuration n with respect to the F conformation is given by AG n (/), the free 
energy difference at a given force / between the duplex containing N — n closed bps and the 
completely closed configuration (F state), 

AG n (/) = AG n (0) + AG^ sRNA (f) + AG*>(/). (1) 

In eq. [l]AG n (0) is the free energy difference at zero force between a hairpin in the partially 
unzipped configuration n and a hairpin in the completely closed configuration; AG^ sRNA (f) is 
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equal to the reversible work needed to stretch the ssRNA strands of the hairpin in configuration 
n (2n opened bases) from a random coiled state to a force-dependent end-to-end distance x a (f); 
and A.G~°(f) is the contribution related to hairpin stem orientation [SHI [57]. An estimation of 
AG n (0) at 1M [Mon + ] can be obtained by using the nearest-neighbor (NN) energy parameters 
widely employed to predict the stability of RNA secondary structures [T5J [TU] . It is given by the 
sum of the stacking contributions of the duplex region, containing N — n bps. The elastic term 
AG s a sRNA (f) is given by [56] 
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The molecular extension of ssRNA, x^ sRNA (f), can be estimated using polymer theory (see 
section 2.4). Finally, the last term in eq. [l] AG*(/), is equal to the free energy of orientation 
of a monomer of length d along the force axis [57J : 



AG*>(/) = fc B Tlog 
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where / is the applied force, k B is the Boltzmann constant, T is the bath temperature, and 
do is the diameter of a double stranded chain, taken equal to 2 nm. 



2.4 Elastic models of ssRNA 

To model the elastic response of ssRNA we employed both the interpolation formula for the 
inextensible Worm Like Chain (WLC) model and the Freely Jointed Chain (FJC) model, which 
give the equilibrium end-to-end distance a; of a polymer of contour length l n stretched at a given 
force /. These models have been mainly tested for long polymers. However, several studies 
indicate that they are generally applicable when the contour length is larger than the persistence 
length. The inextensible WLC is given by: 
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where k B is the Boltzmann constant, T is the bath temperature and P is the persistence 
length [551 [53]. The FJC model is given by 
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(5) 



k B Tj fb 
where b is the Kuhn length. 

There are other models, such as the Thick Chain, that are more general than the WLC or the 
FJC and that have been used to fit the elastic response of biopolymers. Despite of their greater 
complexity, we do not expect a qualitative improvement of our results by using them. 
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2.5 Kinetic Analysis 

We applied Kramers rate theory [BU] to study the kinetics of the transition between states F 
and U. The framework for understanding the effect of an external force on rupture rates was 
first introduced in [6T] and extended to the case where the loading force increases with time 
|62L I63| . The assumption that the transition state does not move under an applied force / can 
be relieved by considering that the effective barrier that must be crossed by a Brownian particle 
is force-dependent, B c g(f). The unfolding and folding rates can be obtained as the first passage 
rates over the effective barrier, 

ku (f) = k exp(-^^, (6a) 
t m i ( B cS (f) - AG N (f) \ 

k F (/) = fco exp I — I . (6b) 

In eq. |6j F was selected as the reference state and AGn(/) has been defined in eq. [I] fco is 
the attempt rate for activated kinetics. The effective barrier B e g(f) can be obtained analytically 
from Kramers rate theory (KT) [BUGIS] (detailed derivation provided in the Supporting Material, 
section S3) as 



B^ T (/) = fc B riog 



(7) 



with h(n) = exp (— ^"'j^ ) • Importantly, the location of the barrier along the reaction 

n'— 

coordinate can be obtained from the first derivatives of B c g(f) with respect to force, 

df 
df 

where x^ s (f) and x^ s (f) are the distances from the effective barrier to the F and U states, 
respectively. The force-dependent fragility parameter fi(f) [36 , 



xf ff (/)-^ ff (/) 

< ff (/)+< ff (/) 

lies in the range [-1:1] and is a measure of the compliance of a molecule under the effect of 
tension. Compliant structures deform considerably before the transition event and are character- 
ized by positive values of fi(f), i.e. x^sif) > x cs(f)- ^ n contrast, brittle structures are defined 
by negative values of At(/), x^{f) < x^ s (f). A given sequence can display different fragilities at 
different force regimes, due to changes in the location of the transition state (TS) with force. 

From the measured transition rates (see section 2.6) we can get estimators for the effective 
barrier B^ F ^ (/) for unfolding and folding using the expressions in eq. JfjJ 
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-^L±L± = -logku(f)+logko, (10a) 
k B T 

B f^i = - log k F ( f) + log kQ + A ^m. dob) 

k B T k B T 
By comparing the experimental estimators of the kinetic barrier 

B is /F) (/) with the effective 

barrier B^ r (f) as predicted by Kramers rate theory (eq. [7]) we can extract the free energy of 
formation of the hairpin AGn(O), the attempt rate ko and the parameters that characterize the 
elastic response of the ssRNA [35J . While k always can be determined by doing this comparison, 
there is a trade-off between the contributions of the elastic response of the ssRNA and the free 
energy of formation of the hairpin. Although this is not strictly true (the stretching contribution 
term is force dependent whereas the free energy of formation term is not) it holds to a very good 
degree. Therefore, if only the free energy of formation of the hairpin is known a priori, then we 
can extract the elastic properties of the ssRNA by matching eqs. [I0]a and [Hj. with eq. [7] On 



the contrary, if we only know the elastic properties of the ssRNA, then we can extract the free 
energy of formation of the hairpin (see Supporting Material, section S4). 

2.6 Data Analysis 

The molecular transitions during unfolding and folding can be identified as force rips in a force- 



distance curve (FDC) [56] . In order to extract the unfolding and folding rates (eq. 6a and 6b) 
from experiments we have collected the first rupture forces associated with the unfolding and 
folding parts of each pulling cycle (Fig. 2 A and 2B). By plotting the number of trajectories 
in which the molecule remained at the initial configuration (F state during the stretching part 
and U state in the releasing part of the cycle) as a function of force N(fi), divided by the 
total number of trajectories Nq, we obtained experimental estimators for survival probabilities 
Pu/Fifi) — N(fi)/No of the U and F states. Moreover, we obtained an experimental estimator for 
the probability densities Pu/F^f) of unfolding and folding first rupture forces by doing normalized 
histograms of both datasets {pu/f(D — AN/(Af x N ), where AN is the number of events in 
the range between / and / + A/). The survival probabilities are related to pjj(F){f) by the 
following equations, 

Pu(f) = 1 - f PuU'W, (Ha) 



J fmin 

P F {f) = l- I 1 ™* PF(f)df . (lib) 
Jf 

If we assume a two-state transition, the time-evolution of the survival probabilities is described 
by the following master equations [55] : 
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dPui / t {t)) = -ku(f(t))Pu(f(t)), (12a) 

dPFi / t (t)) = -k F (f(t))P F (f(t)). (12b) 

With this assumption and the experimental estimators for survival probabilities and densities, 
it is possible to extract the transition rates fc[/(F)(/) from rupture force measurements using 
k u(F)U) = -rpu(F)(f)/Pu(F)(f), with r the pulling speed [13 [53 [Hi] . 

2.7 Salt corrections 

It is interesting to experimentally measure the effect of salt on the free energy of formation 
of nucleic acid hairpins. However, UV absorbance experiments cannot be carried out for this 
particular sequence because its melting temperature is too high to obtain reliable results (see 
Supporting Material, section S5). Therefore, as mentioned in section 2.3, the estimation of 
the free energy of formation of the RNA hairpin at 1M [Mon + ] is obtained using the NN energy 
parameters proposed by [T^JES]- To introduce the effect of monovalent salt concentration [Mon + ] 
we assume a sequence- independent correction gi([Mon + ]) for the free energy of formation of one 
base pair. As the free energy is measured relative to the F state we get for the free energy 
correction of a hairpin with n unzipped bps: 

AG[ Mo,1+l (0) = AGi M (0) - n 5l ([Mon+]) (13) 

where AG„ M (0) corresponds to the free energy of formation of the n th configuration at 1000 
mM [Mon + ] at zero force. In the case of mixed monovalent/Mg 2+ conditions we add a second 
sequence-independent correction term ^([Mon" 1 "]) that captures the effect of Mg 2+ ions on the 
hairpin free energy of formation: 

AG [Mo„+],[M g 2+] (0) = AG 1M (0) _ n5l([M on+]) - n 52 ([Mg 2 +]). (14) 

In what follows, unless stated otherwise, all monovalent and divalent salt concentrations 
[Mon + ], [Mg 2+ ] are expressed in mM units. 



3 Results 

3.1 Effect of force on thermodynamics and kinetics of an RNA hairpin 

We pulled the RNA hairpin at loading rates of 1.8 pN/s and 12.5 pN/s in buffers containing 
different ionic conditions (see Materials and Methods, section 2.2). From the unfolding and 
folding FDC we measured the first rupture forces along many cycles. The resulting probability 
distributions at each pulling speed and ionic condition can be found in the Supporting Material 
(section S6). As previously observed |36j . this hairpin displays a two-state behavior, with force 
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jumps signaling the transition between F and U states (Fig. 2A) and with no evidence of fraying 
or intermediate states (see Supporting Material, section S2). 

The order of magnitude of the resulting rupture forces and hysteresis effects are compatible 
with previous force-melting experiments carried out for other simple nucleic acid structures, like 
P5ab and TAR RNA hairpins [311 131123] or short DNA hairpins [55]. Moreover, results at 1 M 
[Na+] were in significant agreement with solution predictions [TSJ [73]. As expected, hysteresis 
effects strongly depend on the loading rate, being lower at 1.8 pN/s and higher at 12.5 pN/s. 
This can be seen from the experimental distributions of unfolding and folding first rupture forces 
(Fig. 2B), which are closer for pulling cycles performed at 1.8 pN/s. 

3.2 Experiments at different monovalent cation concentrations 

We performed pulling experiments at four different NaCl concentrations (see Materials and Meth- 
ods, section 2.2). We find that the RNA duplex stability increases at higher [Mon + ] concentra- 
tions. For instance, rupture force distributions are displaced to higher forces (Fig. 3A) as we 
increase the concentration of NaCl. The greater duplex stability at higher salt concentrations 
can also be observed as an increase in the mean rupture force with the logarithm of the salt 
concentration (Fig. 3B). The standard deviation of rupture unfolding (folding) forces, that are 
known to be proportional to kgT '/x^g (f) [56] . remain almost constant along the salt range 
explored. That might denote that the position of the TS mediating the unfolding and the folding 
transitions does not depend on [Mon + ], despite the fact that both transitions occur at higher 
forces. In Fig. 3C we see that the unfolding (folding) kinetic rates decrease (increase) with the 
salt concentration, which again shows the stabilizing effect of salt on the RNA hairpin. 



3.3 Experiments at 1M NaCl 

From the current set of NN energy parameters for RNA secondary structures obtained at 1000 
mM NaCl concentration we can predict the free energy of formation AGn(O) for the RNA hairpin 
at this particular condition using the mfold server PH [HI GUI EE] • We get AG^ fold (0) = 63.0 
k^T (although our experiments were performed at 1050 mM of monovalent salt, we do not expect 
significant differences by comparing with the prediction at 1000 mM). By applying the kinetic 
method introduced in section 2.5 we can evaluate the kinetic barrier associated to the unfolding 
reaction B* T {f) and B^ /F) (f) for a given elastic model for the ssRNA and find the one for 
which the theoretical prediction by Kramers rate theory (eq. [7| best matches the experimental 
results (cq. 10). The procedure is shown in Fig. 4A and explained in the Supporting Material 
(section S4). We found the best fit to our data using the inextensible WLC model with persistence 
length P = 0.75 ± 0.05 nm and interphosphate distance a = 0.665 nm/base. The free energy 
of formation obtained is equal to AGn(0) = 65.3 ± 0.3 k^T, in reasonable agreement with the 
aforementioned value for AGN fold (0). 
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3.4 Selection of elastic parameters: experiments with monovalent salts 

In order to know AGn(O) at monovalent ionic conditions different from 1050 mM we need to 
know the effect of salt on the elastic contribution of ssRNA strands. The elastic behavior of 
single-stranded RNA (poly-U) has been studied in single-molecule stretching and fluorescence 
experiments carried out at various [Na + ] concentrations [5j2 EH IZI]- Despite the extremely 
different contour length of the molecules under consideration (we are dealing with 44 bases-long 
chain in contrast to a polynucleotide of 1500-4000 bases in [70]), we take the values for the 
persistence length P proposed in previous stretching experiments and add our value obtained at 
1050 mM. We assume that the elastic properties of ssRNA strands are independent of sequence, 
which can lead to a small error in the values obtained for the elastic properties of ssRNA strands 
in the case of sequence-dependent behavior. In fact, a sequence-dependent elastic behavior for 
ssDNA strands was previously considered as a possible explanation for the specific salt corrections 
found for the NN energy values obtained from unzipping experiments [37] . In the inset of Fig. 
5 we plot the persistence lengths versus the monovalent salt concentration. To fit the data, we 
employ the following dependence of P on the Debye screening length \ V D: 

P~\ u n oc (15) 

D [Mon + ]«V 2 V ' 

No unique scaling law can be derived, as the value of v can be 1, 2 or <1, depending on the 
polymer properties [75]. We can fit all data using v = 1 but also letting v be a free parameter, 
which results in a value of v = 0.6 ± 0.06. We then interpolate both fits, in order to infer the 
values of P in our experimental conditions of [Mon + ] concentration (50, 150, 550 and 1050 mM), 
as shown in Table 1. 

3.5 Monovalent salt correction to the free energies of formation of the 
RNA hairpin 

Having obtained the elastic parameters that allow us to appropriately describe the elastic re- 
sponse of ssRNA strands at different [Mon + ] concentrations, we still need to characterize the 
effect of salt on the energies of formation of the RNA hairpin at each intermediate configuration 
n. It is generally assumed a sequence-independent correction to the free energies of formation of 
nucleic acids duplexes [73, 74, 75 . However, we have previously shown that a sequence-dependent 
salt correction to the NN energy parameters of DNA improves the free energy prediction of both 
unzipping and melting experiments [37] . Related to this, it has been found that cation con- 
centration affects RNA stability in a sequence-dependent manner [75] . In the absence of RNA 
sequence-specific parameters available, we adopted a sequence-independent salt correction (eq. 



13) given by gi([Mon + [) = mlog([Mon + ]/1000), where [Mon + ] is expressed in mM units. As we 
will see, there are experimental and theoretical evidences that support the logarithmic effect of 
monovalent ions to the stability of nucleic acid hairpins. 
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Using this correction, the variation of AG n (0) with monovalent salt concentration depends 
strictly on the value of the constant m. In order to derive m from our data, we compared 
the estimators of -B e ff(/) obtained experimentally (B^\f) and B^P(f) in cq. 10) with the 
theoretical prediction {B^f (f) in eq. [7]) at different values of m. In Fig. 4A-D, we see the 
correspondence between theory and experiments at each monovalent ion concentration. For all 
salt concentrations, we found the best agreement at m — 0.10±0.01 kcal/mol. This value agrees 
with the sequence-independent salt correction reported for DNA duplex oligomers in melting 
experiments, m = 0.110 ± 0.033 kcal/mol [731 177j . and in unzipping experiments of polymeric 
DNA, m = 0.104 ± 0.010 kcal/mol [37]. Fig. 3D summarizes all the results. At a given force we 
see that the height of the kinetic barrier increases with salt concentration, which again indicates 
that salt increases kinetically the stability of the RNA structure. 

In Fig. 5 we show the dependence of the measured AGn(O) of the RNA hairpin on the 
monovalent ion concentration. As expected from earlier observations on DNA [351 ESI US] and 
from the application of counterion condensation theory to interpret polyelectrolyte effects on 
equilibrium involving highly charged, locally rod-like polyelectrolytes [40l [80l HH [82] , we observe 
an approximately linear dependence of RNA duplex stability on the logarithm of monovalent 
salt concentration. Interestingly, our data can also be well-described by the empirical expressions 
derived in [43 , 44 , 45j , where the TBI model is used to predict the hairpin free energies at different 
ionic conditions (see Supporting Material, section S7). 

By deriving the effective barrier as a function of force we can measure the distance of the 
TS to the F and U states, x^ s (f) and x^ s (f) (eqs. 8a and 8b), and the fragility fi(f) of the 



molecule as a function of the applied force (eq. [9|. Fig. 6 shows the two extreme cases with 50 
and 1050 mM [Mon + ] (continuous and dashed lines respectively). In panel A we observe that 
the location of the TS changes as a function of force. The same trend is observed for the fragility 
in panel B, where the experimentally measured points, the predicted force-dependent fragility 
(black curves), and the expected values of the fragility for all possible locations n of the TS along 
the stem of the hairpin are represented (horizontal grid, right scale). At low forces, the TS is 
located near the loop, n = 19 ±1 (Fig. 6A, dark gray curves). At intermediate forces that depend 
on salt concentration the TS moves to the stem region n = 6± 1 (Fig. 6A, gray curves). At large 
forces the TS has disappeared. These results are in agreement with previous findings using the 
same hairpin sequence [3J>]. Moreover, we see that at higher monovalent salt concentrations the 
locations of the different TS mediating unfolding and refolding are the same (n ~ 18, 6 or for 
low, intermediate and high forces respectively) but shifted to larger forces. These results agree 
with the Hammond's postulate [83]: at increasing [Mon + ] the F state is increasingly stabilized 
while the TS is shifted towards the U state; simultaneously, as force increases the TS approaches 
the F state. 
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3.6 Experiments in mixed monovalent/Mg 2+ conditions 

We have also performed pulling experiments in mixed monovalent /Mg 2+ buffers, containing a 
fixed concentration of Tris + ions (50 mM) and varying concentrations of Mg 2+ (see Materials 
and Methods, section 2.2). The rupture force distributions for all mixed monovalent/Mg 2+ 
conditions can be found as Supporting Material (section S6). We found two regimes in the 
behavior of the average rupture forces for unfolding and folding processes along the range of 
[Mg 2+ ] experimentally explored. Below 0.1 mM [Mg 2+ ], there is no significant difference between 
control (no Mg 2+ added) and magnesium-containing conditions (Figs. 7A and 7B). However, 
at higher magnesium concentrations, we found a linear dependence of average rupture forces 
with the logarithm of [Mg 2+ ] (Fig. 7A and 7B). Interestingly, Owczarzy et al. have made a 
similar observation in DNA melting experiments done in mixed monovalent/Mg 2+ conditions 



[5"T] . They found that the ratio R = J [Mg 2+ ]/[Mon + ] (both salts in molar units) is a convenient 
parameter to determine whether divalent or monovalent ions are dominant on duplex stability. 
If R is less than 0.22, then monovalent ions are dominant and the presence of Mg 2+ can be 
ignored. In our experiments, R = 0.0632 and R = 0.2 for the 0.01 mM and 0.1 mM [Mg 2+ ] 
conditions, respectively. As in the case of pure monovalent ion conditions, the standard deviation 
of rupture forces remains almost constant and we also observed a linear dependence of log &[/(/) 
and logfcp (/) on the applied force for the different [Mg 2+ ] tested (Fig. 7C). 

In order to obtain the free energy of formation AGn(0) of the RNA hairpin at different 
magnesium concentrations, we employed the empirical expression derived in [441 45] by applying 
the TBI model to predict the RNA helix stability in mixed monovalent /Mg 2+ ionic conditions (see 
Fig. 9, inset). Using this mixed salt correction, it is possible to obtain the sequence-independent 
correction of one base pair <? 2 ([Mg 2+ ]) at any mixed salt condition using eq. 
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32 ([Mg 2+ ]) = i (AG^ fold (0) - AG™(0) - 7Vmlog([Mon+/1000])) . (16) 



From this expression, we can extract the value of AG l ™ on ] ' [Mg ] (0) in eq. Il4l for any 
intermediate state n. By varying the mixed salt-dependent values for the persistence length P of 
the ssRNA for each [Mg 2+ ], we can now determine the value of P that results in better agreement 
between the predicted effective barrier B^ r (f) and our experimental estimations (Fig. 8 and 
Supporting Material, section S4). All the results are summarized in Fig. 7D, where we can see 
that the stability of the hairpin increases with magnesium concentration. The dependence of P 
on [Mg 2+ ] is shown in Fig. 9. Table 2 summarizes the results obtained for the persistence length 
P and the attempt frequency ko. The position of the TS varied with [Mg 2+ ] in a way similar 
to what we found for [Mon + ]. In a specified force the TS mediating the unfolding and folding 
transitions is shifted toward the U state as the [Mg 2+ ] is raised (Fig. 10A), in agreement with 
the Hammond's postulate jH3|- The force-dependence of the position of the TS with respect to 
the F state x^ s (f) and the hairpin fragility //(/) are similar in both monovalent and mixed ionic 
conditions (Figs. 6A, B and 10A, B). At low forces the TS is located near the loop, whereas at 
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intermediate forces it is located in the stem region n — 6 ± 1 (Fig. 10A, B). 

4 Discussion 

The effect of monovalent ion concentration on DNA stability has been extensively studied and 
there is a variety of empirical salt corrections available in the literature (HH [85] . There is no gen- 
eral agreement about the accuracy and limitations of use of salt corrections in terms of sequence 
length and range of salt concentrations [53]. Recently, we have reported 10 NN salt correction 
parameters for prediction of DNA duplex stability derived from single-molecule experiments [37J . 
However, there is no equivalent study on RNA duplexes and the experimental data available for 
the salt effects on RNA duplex stability are limited to short sequences that display a two-state 
behavior [75] . Different polyelectrolyte theories try to characterize the interaction between coun- 
terions and nucleic acids to study ionic effects to the molecular stability. The most accepted 
are mean field theories such as the Poison-Boltzmann and the counterion condensation theories 
[8"Ul 15T] . Recently, the Tightly Bound Ion (TBI) model has been proposed [33J . It incorporates 
correlation and fluctuation effects for bound ions, and has been extended to treat RNA helices 
under mixed monovalent /divalent salt conditions |44l 145) . It was shown that the TBI improves 
the prediction of the stabilities of RNA duplexes smaller than 15 bp [44l 186] , 

Here, we performed a detailed characterization of the effect of monovalent and mixed mono- 
valent /magnesium concentrations on the stability of a RNA hairpin containing a stem of 20 
bp by mechanically unfolding and folding the molecule using optical tweezers. The results we 
have obtained can be very well described by the empirical formulas derived from the TBI model 
for predicting the stabilities of RNA hairpins in monovalent ion and mixed monovalent/Mg 2+ 
buffer conditions. A comparison of the counterion condensation theory 87J is provided in the 
Supporting Material (section S8), and although an acceptable accuracy is found at monovalent 
salt conditions, we observe that correlation effects become important in the presence of Mg 2+ 
ions and therefore the TBI gives improved predictions. 

When a hairpin is melted mechanically the U state is subjected to a tension. In standard 
optical melting and calorimetric experiments, the U state is under zero force. Therefore, to ap- 
propriately compare both phenomena, it is necessary to discount the energy required to stretch 
the single stranded unfolded state in mechanical experiments. However, we have limited infor- 
mation about the elastic properties of single stranded RNA as compared to DNA. Despite the 
chemical similarity between both polymers, the RNA-specific 2'-OH group favors a conformation- 
ally restricted C3' — endo structure [55]- Moreover, in previous single- molecule studies [5§1 ITT] , 
the elastic response of a single-stranded poly-U RNA was better described by a WLC model 
rather than by a FJC model. Based on these considerations, we decided to use the numbers 
derived from single-molecule stretching experiments for the elastic behavior of ssRNA molecules 
(as explained in section 3.2) to appropriately account for the energetic contribution of stretching 
the ssRNA strands during mechanical unfolding. 



12 



By applying the Kramers rate theory and correctly accounting for the elastic contribution of 
ssRNA stretching we were able to obtain the free energy of formation of the RNA hairpin from 
pulling experiments done in non-equilibrium conditions. This allowed us to obtain a large number 
of trajectories at different ionic conditions in a feasible timescale. In this way, we obtained the 
effective barrier of the unfolding reaction as a function of force for different ionic conditions. 

The ability of Mg 2+ ions to stabilize RNA structures at much lower concentrations than 
monovalent ions was recognized almost forty years ago [89]. In fact, by plotting the values of 
the free energies of formation versus salt concentration (expressed in mM units) we can collapse 
data for both types of salt into a single master curve by multiplying [Mg 2+ ] by a factor 100 (Fig. 
11A). This effect has been previously observed [231 ESJ |3H [501 IE] and can be explained using 
the counterion condensation theory which account for strong correlations between counterions 
and polyelectrolytes [33 [HU E2 E2] ■ A similarity (but not a data collapse) is obtained for the 
persistence length values of ssRNA (Fig. 11B). It can be interpreted as the screening effect of the 
counterions present in solution: the longer the screening length, the stiffer the molecule. This 
equivalence found between the non-specific binding of monovalent and divalent salts with nucleic 
acids might be useful to develop new biochemical assays in situations where divalent ion-specific 
interactions could be detrimental |93) . 

As a further step, it would be very interesting to extend this study by mechanically un- 
folding an RNA hairpin containing a Mg 2+ -specific binding site using our experimental setup. 
Eventually, the energetic contributions of different specific binding sites could be dissected and 
incorporated into current models of RNA structure prediction. 

5 Supporting material 

Supporting material is available at NAR Online. 
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7 Figure Legends 
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Figure 1. Experimental setup and RNA sequence. (A) A single RNA hairpin is 
attached to two polystyrene microspheres through RNA/DNA heteroduplexes used as handles. 
The anti-digoxigenin antibody-coated microsphere is optically trapped while the streptavidin- 
coated microsphere is positioned at the tip of a micropipette by air suction. (B) RNA hairpin 
sequence. 
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Figure 2. Pulling experiments. (A) A few pulling cycles for the RNA hairpin showing the 
unfolding (red) and folding (blue) trajectories. (B) Experimental distribution for the unfolding 
first rupture forces at 1.8 pN/s (red full squares) and 12.5 pN/s (blue full circles), and for the 
folding first rupture forces at 1.8 pN/s (red empty squares) and 12.5 pN/s (blue empty circles). 
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Figure 3. Kinetic analysis of experiments at varying [Mon + ]. (A) Experimental 
distribution of the unfolding rupture forces in buffers containing 50 mM (red), 150 mM (green), 
550 mM (blue), and 1050 mM (magenta) [Mon + ]. These experiments were done at a loading rate 
of 1.8 pN/s. (B) Average rupture forces (fu(F)) an( A standard deviations cr[/(i?) as a function of 
monovalent cation concentration at loading rates of 1.8 pN/s (red) and 12.5 pN/s (blue). Full 
symbols refer to unfolding and empty symbols to folding. (C) Log-Linear plot of the transition 
rates versus force. Experiments were done at 50 mM [Mon + ] for loading rates of 1.8 pN/s (dark 
red) and 12.5 pN/s (red), at 150 mM [Mon + ] for loading rates of 1.8 pN/s (dark green) and 12.5 
pN/s (green), at 550 mM [Mon + ] for loading rates of 1.8 pN/s (dark blue) and 12.5 pN/s (blue), 
and at 1050 mM [Mon + ] for loading rates of 1.8 pN/s (dark violet) and 12.5 pN/s (magenta). 
(D) Dependence of the effective barrier B c s(f) on force at different [Mon + ]. Color code as in 
(C). 
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Figure 4. Determination of the salt correction parameter m. Experimental evalua- 
tion of the sequence- independent correction parameter m from gi([Mon + ]) = mlog([Mon + ]/1000) 
using eq. 
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Estimators of B^g (/) obtained experimentally (eqs. 10 1 were compared with 



the expected B^f(f) (eq. [7]) profiles for different values of m (0, 0.1, 0.2, 0.3, 0.4 and 0.5 from 
top to bottom). Red (green) points are the experimental estimators B^\f) at a pulling rate of 
12.5 (1.8) pN/s. Blue (magenta) points are the experimental estimators of B^P(f) at a pulling 
rate of 12.5 (1.8) pN/s. Light blue lines are the B^ r (f) profiles for values of m not matched, and 
black lines are the experimental estimators of B c tf(f) that match with experiments. Application 
of the method to experiments done at 1050 mM [Mon + ] (A), 550 mM [Mon + ] (B), 150 mM 
[Mon+] (C), and 50 mM [Mon+] (D). 
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Figure 5. Free energy of formation of the RNA hairpin as a function of [Mon + ]. 

Main panel: Free energy obtained experimentally (squares), using the logarithmic dependence 
with salt concentration given by gi([Mon + ]) (dashed line) and using the TBI model (continuous 
line) Inset: Persistence lengths P obtained from the application of the Thick Chain model to 
published experimental data for poly-U RNA stretching in buffers containing 5, 10, 50, 100, 300, 
and 500 mM of [Na + ] (squares) [69j[70]. We have included the value of the persistence length 
P that we obtained in this study at 1050 mM [Mon + ] (empty circle). Two different fits to data 
were done from eq. 15 a fixed value v = 1 (red) and v as free parameter v — 0.60 ± 0.06 (blue). 
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Figure 6. Barrier location and mechanical fragility at 50 mM and 1050 mM 
[Mon + ]. (A) Force-dependence of the barrier position measured with respect to the F state, 
x cs(f)- Continuous gray line is the WLC prediction of the molecular extension when n = 
19 or n — 6 bps are unzipped at 50 mM [Mon + ], and dashed gray line corresponds to the 
WLC prediction when n = 19 or n = 6 bps are unzipped at 1050 mM [Mon + ]. As seen, at 
an intermediate value of forces n = 6 coincides with the TS for both ionic conditions. (B) 
Dependence of fragility ^(/) with force. Gray lines indicate the value of the fragility for different 
locations n of the TS along the stem. Continuous black lines are the theoretical prediction using 
Kramers rate theory for data at 50 mM [Mon + ], and dashed black lines for data at 1050 mM 
[Mon + ]. Blue and green points are the experimental evaluation of x^ s (f) and for folding 
and unfolding data collected at 50 mM [Mon + ]. Red and purple points are the experimental 
evaluation for folding and unfolding at 1050 mM [Mon + ]. 
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Figure 7. Kinetic analysis of experiments at varying [Mg 2+ ]. (A) Experimental 
distribution of the unfolding rupture forces in buffers containing 0.00 mM (red), 0.01 mM (green), 
0.1 mM (blue), 0.5 mM (magenta), 1 mM (cyan), 4 mM (orange), and 10 mM (black) of [Mg 2+ ] 
and 50 mM of monovalent cations. These experiments were done at a loading rate of 1.8 pN/s. 
(B) Average rupture forces (fu(F)) and standard deviations <Ju(f) obtained in experiments done 
at different [Mg 2+ ] and at loading rates of 1.8 pN/s (red) and 12.5 pN/s (blue). Full symbols 
refer to unfolding and empty symbols to folding. (C) Log-Linear plot of the transition rates 
versus force. Experiments were done at 0.00 mM [Mg 2+ ] for loading rates of 1.8 pN/s (dark red) 
and 12.5 pN/s (red), at 0.01 mM [Mg 2+ ] for loading rates of 1.8 pN/s (dark green) and 12.5 
pN/s (green), at 0.1 mM [Mg 2+ ] for loading rates of 1.8 pN/s (dark blue) and 12.5 pN/s (blue), 
at 0.5 mM [Mg 2+ ] for loading rates of 1.8 pN/s (dark violet) and 12.5 pN/s (magenta), at 1 
mM [Mg 2+ ] for loading rates of 1.8 pN/s (dark cyan) and 12.5 pN/s (cyan), at 4 mM [Mg 2+ ] 
for loading rates of 1.8 pN/s (dark orange) and 12.5 pN/s (orange), and at 10 mM [Mg 2+ ] for 
loading rates of 1.8 pN/s (gray) and 12.5 pN/s (black). (D) Dependence of the effective barrier 
B fi(f) on force at different [Mg 2+ ]. Color code as in (C). 
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Figure 8. Determination of the persistence length of ssRNA at varying [Mg 2+ ]. 

Estimators of B^ F \f) obtained experimentally were compared with the expected B^{f) 



profiles for different values of P (0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4 and 1.5 nm from top 
to bottom) using eq. [7] and eqs. 10 Red (green) points are the experimental estimators B^P(f) 



at a pulling rate of 12.5 (1.8) pN/s. Blue (magenta) points are the experimental estimators of 
B^g '(/) at a pulling rate of 12.5 (1.8) pN/s. Light blue lines are the B^f{f) profiles for values 
of m not matched, and black lines are the experimental estimators of B c g(f) that match the 
experiments. Application of the method for experiments done at 0.01 mM [Mg 2+ ] (A), 0.1 mM 
[Mg 2 +] (B), 0.5 mM [Mg 2 +] (C), 1 mM [Mg 2+ ] (D), 4 mM [Mg 2+ ] (E), and 10 mM [Mg 2 +] (F). 
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Figure 9. Dependence of the persistence length P on [Mg 2+ ]. Main panel: Experi- 
mental persistence length versus [Mg 2+ ]. Inset: Dependence of the free energy of formation of 
the RNA hairpin AG^ BI (0) on [Mg 2+ ] with fixed 50 mM [Mon+] obtained using the TBI model 
|44l|45]. Black points are the values of the free energy of formation that we used for our analysis. 
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Figure 10. Barrier location and mechanical fragility at 0.01 mM and 10 mM 

[Mg 2+ ]. (A) Force-dependence of the barrier position measured with respect to the F state, 
x cs(f)- Continuous gray line is the WLC prediction of the molecular extension when n = 19 
or n = 6 bps are unzipped at 0.01 mM [Mg 2+ ], and dashed gray line corresponds to the WLC 
prediction when n = 19 or n = 6 bps are unzipped at 10 mM [Mg 2+ ]. At an intermediate range of 
force the TS coincides with n — 6 for both ionic conditions. (B) Dependence of fragility n(f) at 
0.01 mM and 10 mM [Mg 2+ ]. Gray lines indicate the value of the fragility for different locations 
n of the TS along the stem. Continuous black lines are the theoretical prediction using Kramers 
rate theory for data at 0.01 mM [Mg 2+ ], and dashed black lines for data at 10 mM [Mg 2+ ]. Blue 
and green points are the experimental evaluation of x^ s (f) and n(f) for folding and unfolding 
data collected at 0.01 mM [Mg 2+ ]. Red and purple points are the experimental evaluation for 
folding and unfolding at 10 mM [Mg 2+ ]. 
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Figure 11. Comparison between [Mon + ] and [Mg 2+ ] results. (A) Free energy of 
formation of the RNA hairpin at different salt conditions. Magnesium concentrations have been 
multiplied by 100 along the horizontal axis. (B) Persistence length values for the ssRNA hairpin 
at different salt conditions. Magnesium concentrations have been multiplied by 100 along the 
horizontal axis. 
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8 Tables 



[Mon+] (mM) 


P (nm) 


logfc (1/s) 


AG N (0) (k B T) 


1050 


0.75±0.05 


10.9±0.4 


65.1±0.3 


550 


0.82±0.02 


10.5±0.4 


64.0±0.4 


150 


1.01±0.01 


11.3±0.5 


59.4±0.4 


50 


1.27±0.03 


12.4±0.4 


54.0±0.4 



Table 1: Parameters obtained from experiments at different [Mon + ]. ssRNA persis- 
tence length P, logfco, and free energy of formation (AGn(0)) for the RNA hairpin at different 
monovalent ion concentrations. 



[Mg 2 +] (mM) 


P (nm) 


logfc (1/s) 


AG^(0) (k B T) 


g 2 (kcal/mol) 


0.00 


1.27±0.03 


12.40±0.40 


55.58 


0.000±0.005 


0.01 


1.50±0.15 


12.05±0.30 


55.50 


-0.007±0.005 


0.10 


1.25±0.10 


11.45±0.30 


55.55 


-0.005±0.005 


0.50 


0.90±0.15 


11.50±0.30 


57.06 


0.0393±0.005 


1.00 


0.80±0.10 


11.40±0.50 


58.63 


0.0858±0.005 


4.00 


0.75±0.10 


11.15±0.50 


62.60 


0.2033±0.005 


10.0 


0.75±0.10 


10.40±0.50 


64.77 


0.2678±0.005 



Table 2: Parameters obtained for experiments at different [Mg 2+ ]. Persistence length 
for ssRNA P, logfc , theoretical predictions for the free energies of formation based on the 
TBI model, AG^ B/ (0), and sequence-independent correction g 2 ([Mg 2+ ]) for the RNA hairpin 
at different magnesium concentrations. A fixed concentration of 50 mM [Mon + ] was used in all 
ionic conditions. 
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Supporting Material 
9 Molecules studied 



[Mon+] (mM) 


[Mg 2 +] (mM) 


loading rate (pN/s) 


Molecules 


Total cycles 


50 





1.8 


4 


95 


50 





12.5 


2 


376 


150 





1.8 


6 


292 


150 





12.5 


4 


329 


550 





1.8 


5 


163 


550 





12.5 


3 


501 


1050 





1.8 


9 


185 


1050 





12.5 


9 


405 


50 


0.01 


1.8 


5 


146 


50 


0.01 


12.5 


5 


386 


50 


0.1 


1.8 


7 


374 


50 


0.1 


12.5 


9 


1434 


50 


0.5 


1.8 


2 


112 


50 


0.5 


12.5 


2 


533 


50 


1.0 


1.8 


4 


205 


50 


1.0 


12.5 


6 


2183 


50 


4.0 


1.8 


7 


385 


50 


4.0 


12.5 


7 


1112 


50 


10.0 


1.8 


7 


190 


50 


10.0 


12.5 


3 


1189 



Table 3: Number of molecules and total cycles measured at each ionic salt condition. 
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10 Study of fraying 



The phenomenon of "fraying" at the ends of DNA and RNA duplexes can potentially interfere 
with both solution and single-molecule measurements of DNA and RNA stabilities, and it was 
previously suggested that its effects should be introduced in data analysis [55] . In order to check 
if fraying has an important role in our sequence, we computed the released (absorbed) molecular 
extension Ax m in the unfolding (folding) process. This can be done using the expression: 

Ax m = (17) 

Keff 

where A/ is the force jump measured along the force-distance curve (FDC) and fc e ff is its slope 
before the transition. The change in the molecular extension is also equal to: 

Ax m = x N (f) - x n (f) (18) 

where x^{f) is the equilibrium end-to-end distance of the unzipped hairpin evaluated at the 
unfolding/folding force (N — 20 in this case); x n (f) is the projection of the folded hairpin along 
the force axis; and n denotes the number of open/frayed base pairs in the F state. Ideally, in the 
absence of fraying, n = 0. However, in presence of fraying we should find the value n > such 



that cq. (17) and (18 1 give the same change in molecular extension. 



From the FDC we obtain Af = 1.2 ± 0.1 pN and fc off = 0.0625 ± 0.0146 pN/nm which, using 



eq. (171, gives Ax m = 19 ± 2 nm for any salt concentration. On the other hand, we evaluate 



xn(J) and x„ = o(/) using the elastic properties summarized in Tables 1 and 2 at the measured 



unfolding/folding forces. Using eq. (18) for n = 0we obtain that predicted values for Ax rn lie in 



the range between 18.2 and 19.7 nm. Therefore, we conclude that fraying is not important for 



the molecule under study because the experimental evaluation of Ax m (eq. ( 17 )) is in agreement 
with the estimation of the released molecular extension for N = 20 and n = 0. 

The effect of fraying has been proved to play an important role in former single-molecule 
stretching experiments, like in Woodside et. al. |55j . To understand when fraying is relevant in 
RNA or DNA hairpins we can take a look to sequences at the beginning of the stem as a higher 
GC-content makes the structure more stable. For instance, the RNA hairpin studied here starts 
with 5'-GCG-3', whereas most sequences studied by Woodside et al. [55J start with 5'-GAG-3' 
(except two sequences that start with 5'-TAT-3' and 5'-AAG-3'). In Table [I] we compute the 
free energy difference AG\(f) at 1 M NaCl between the completely folded conformation (n = 0) 
and the frayed configuration with one open distal base pair (n — 1) at different values of force 
for our RNA hairpin and molecule 20R55/4T in [55]. At zero force, AG^°' R55 ^ 4T (0) is below 3 
k^T and consequently thermal fluctuations can overcome the energetic barrier and the frayed 
conformation can take place. However, in the case of our hairpin AG\ is too high for thermal 
fluctuations to overcome the energetic barrier. The same trend is observed at 10 pN. At 20 pN 
fraying is irrelevant because both molecules are in the unfolded state. 



34 



Force (pN) 


AG 20K55/4T {kBT) 


AGf NA (k B T) 





2.67 


6.45 


10 


1.08 


4.86 


20 


-1.34 


2.45 



Table 4: Free energy differences between frayed (n = 1) and completely closed structures (n = 0) 
of 20R55/4T hairpin [55] and our RNA hairpin at 1 M NaCl and different forces. 



Motivated by solution measurements [531 12Z] , we also considered the possibility of fraying at 
the opposite end of the stem, on the base pairs closest to the loop. However, we found free energy 
differences of more than 10 k-gT between these configurations and the closed configuration, and 
consequently we conclude that these "frayed" configurations are not affecting our results. 

Based on these considerations, we conclude that fraying plays a rather minor role (if any) on 
the thermodynamics and kinetics of folding/unfolding of the RNA hairpin and can be neglected. 



11 Derivation of the effective barrier B^(f) 

Here we derive the analytical expression for the effective barrier of a one-dimensional free energy 
landscape based on the work by Hyeon and Thirumalai |65j . 



V(x) A 




a b x 
Figure 1: Brownian particle in a double well potential. 

Suppose the system sketched in Fig. [TJ where a Brownian particle is subject to the one 
dimensional potential V(x). The time evolution of the probability density function (pdf) p{x,t) 
to find the particle at the position x at time t follows the Fokker-Planck equation [551 E3] : 



dp(x, t) 
dt 



= D 



D 



d 
Ox 
d 



d_ 
Ox 



dx 

£ F p \p(x,t)] 



1 dV(x) 
k^T dx 

l d y w 
dx 



p(x,t) 
p(x,t) 



(19) 



Where D is the diffusion coefficient and Cpp is the Fokker-Plank operator. If we suppose 
that the particle initially is located at x = a, p(x 7 0) = S(x — a), the formal solution of equation 
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((19) is: 

p(x,t) = e tCpp 5{x- a) (20) 

We want to evaluate the average time r it takes to the Brownian particle to jump the kinetic 
barrier located at x — b. In order to simplify the following calculations, we suppose that there 
are absorbing conditions at x — b: once the particle reaches the maximum sketched in Fig. |T| it 
always goes to the right well. The probability to find the particle in the left well of the potential 
V{x) (x G [—00,6]) at time t, also known as the survival probability S(t), can be defined as: 

S(t) = f dx p(x,t) (21) 
J —00 

The time derivative of the survival probability is equal to the time survival pdf p(t), that is, 
the pdf of the time it takes to the Brownian particle to cross the barrier located at x = b. 



S(t + dt) - S(t) = -p(t)dt => p{t) = - (22) 



Therefore, the mean first passage time r can be calculated as: 



poo 

t= dtt p(b,t) 
Jo 



dt 5(6, t) 

/b 
dx p(x, t) 
-00 

/b 
dx e tCrp 6(x-a) 
-OO 

dt / dx 5(x - a)e tCFP \ 

J — 00 

dte tc Frl (23) 

In order to obtain this expression we integrated by parts and used that the adjoint operator 
satisfies f(x)£[g(x)] — g(x)C^ [f(x)]. If we apply the adjoint Fokker-Plank operator £* FP at 
both sides of equation ((23|) we obtain: 



3G 



I -FP T 



dt C FP e tc ^l 



00 , de tc 'FP 
dt — - — 1 



dt 



1 



(24) 



Which gives a differential equation for the survival time r that depends on the adjoint Fokker- 
Plank operator. 

It can be demonstrated that: 



C FP = De«B 



Yisl de k s T d 



dx dx 

and we can write the following differential equation: 



C Fp r(x)=De^ — (e —r{x) 



-1 



(25) 



(26) 



In order to solve eq. ( 26 ) we use the absorbing boundary condition r(6) = and the reflecting 
boundary §^\ x = a = 0. 

1 



t{x) 



D 



dy e fc B T / dxe fcnT 



(27) 



Once we have evaluated the average survival time of the Brownian particle, we want to 
evaluate the effective barrier of the potential V(x). Using the phenomenological Arrhenius 
approach |95j that considers that the survival time depends on the exponential of the barrier we 
can write: 



r ~ e k v T =>• B cS = k B T\og (t/t ) 



B eS = k B T log 



vjy) f y vj*) 
dy e k B T I dxe k B T 



(28) 
(29) 



t is related to the diffusion time of the particle the along x axis. By discretization of equation 



(29) we obtain the expression (7) of the main paper, 



B eS (f) = k B T log 



N 



ri=0 



(30) 



Where AG n is the potential energy V(x) and where we considered that Dtq ~ 0(1). 
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12 Sensitivity of the data analysis 



An experimental estimation of the kinetic barrier is obtained from the measured transition rates 
kjj(f) and k F (f). From the unfolding transition rate, the estimator of the kinetic barrier is given 
by: 



k B T 



log fc - log &[/(/) 



(31) 



Where log fco is a constant (equal to the logarithm of the attempt rate at zero force for the 
activated kinetics) and log ku(f) is estimated from the measured unfolding rupture forces (see 
section 2.6 in the main paper). 

On the other hand, from folding rupture forces the estimation of the kinetic barrier is: 



k B T 



= log fc - log k F (f) + 
= log fco - log k F (f) + 
= log fc - log k F (f) + 



AG N (/) 




k B T 




AG N (0) 


AG ssRNA { 


k B T 


k B T 


AG N (0) 


ff ssRNA 
JO X N 


k B T 


k B T 



(/) , AG*(/) 



(f')df 



k B T 



log 



smh ( 


'fdo\\ 


. fd ' 


yk B Tj 



(32) 



Where log fco is the same constant as in eq. (31 1 and log kp (/) is obtained from the measured 



a salt-dependent persistence length P); and the last term, log 



is evaluated 



folding rupture forces. The term AG^(Q)/k B T is another constant equal to the free energy of 
the RNA hairpin at zero force; Jq x^ RNA (f')df is a force dependent term evaluated according 
to the model used to describe the elastic response of ssRNA (here we use the WLC model with 

using do = 2.0 nm at any salt condition for the given force / (see section 2.3). 

Typically, log fco is unknown and either AGN(0)/k B T is unknown and P is known (here, for 
monovalent salt conditions) or AG at(0) jk B T is known and P is unknown (for mixed monovalent/Mg 2+ 
conditions). 



12.1 Sensitivity of the method at determining AGn(O) 

For a given value of the persistence length P we can evaluate the kinetic barrier from experimental 
unfolding/folding rupture forces by ignoring the unknown constants (log fco and AGAf(0)/fcBT) 
and we obtain the result shown in Fig. [2]A Error bars are evaluated using the bootstrap method. 

In order to determine the constant AGjsi(0)/k B T, i- e. the free energy of formation of 
the RNA hairpin, we impose the continuity of the kinetic barrier in folding and unfolding data 
(Fig. [2j3) . The error committed in the evaluation of AGjv(0)/fcBT mainly depends on the good 
agreement of the overlapping between [B^p (/) — log fco^j and (^B^ (/) — log fco^j and their error 
bars. In the example provided in Fig. [2j3 the best overlapping is found at AGat(O) = 64 k B T. 
In the insets we see that the continuity requirement worsens for values of AGat(O) as close as 
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64.4 or 63.6 k^T. As a consequence, we estimate AGjv(O) = 64.0 ± 0.4 k^T (see Table 1 in the 
main document). 



A B 




/ (pN) / ( P N) 



Figure 2: Determination of AG N (0). Data obtained at 1.8 pN/s and 550 mM [Mon+]. (A) 
Evaluation of B {u \f) - logfc (red squares) and £ (i?) (/) - lo 8' fc o - AGAr(0)/fc B T (blue circles). 
(B) AGjv(0)/fc B r is obtained using the overlapping of the kinetic barrier with force. Once 
unfolding and folding data overlap, the resulting experimental curve is equal to B e g(f) — logfc . 
Insets: The value of AGjv(0)/fcB? 1 can be overestimated (top) or underestimated (bottom) so 
that the continuity requirement is not well-satisfied. 



12.2 Sensitivity of the method at determining logk 

In order to estimate the attempt rate at zero force we need a theoretical model for the kinetic 



barrier. In the case of this work we use the Kramers theory (see section 11). Once we have 
evaluated £? e ff(/) — logfco, logfco is obtained by overlapping the theoretical model to the exper- 
imental results, as shown in Fig. [3] The sensitivity in the determination of log fco is similar to 
the sensitivity in determining AGjy(O): at logfco=10.5 we find the best match, and for values 
0.4 greater or smaller the result significantly worsens (insets in Fig. |3|3). 

12.3 Sensitivity of the method at determining P 



The elastic contribution to B e g(f) only applies to folding data, eq. (32), and regulates the slope 
of the experimental estimation of the kinetic barrier in the range of experimentally measured 
folding forces. In this work we determine the persistence length P of ssRNA by comparing the 
profile of the kinetic barrier estimated from experimental data with the one evaluated using the 
Kramers theory (see Fig. |4| . The sensitivity at determining P is limited by the determination of 
AGjv(0)/fcBT, logfco an d the error bars estimated for Sgjj (/). Therefore, the is a feedback in 
the determination of the AGat(0)//cbT, logfco and P: optimal values are those that give a best 
fit between experimental and theoretical estimations of the kinetic barrier. 
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/ (PN) / ( P N) 

Figure 3: Determination of log k . Data obtained at 1.8 pN/s and 550 mM [Mon + ]. (A) Black- 
straight line is the theoretical evaluation of the kinetic barrier using Kramers theory (section 
S3 and section 2.5 in the main document) and blue circles (red squares) are the experimental 
estimation of the kinetic barrier without the contribution of log fco using folding (unfolding) 
rupture forces. (B) log fco is obtained overlapping the experimental data to the theoretical curve. 
Insets: The value of log fco can be underestimated (top) or overestimated (bottom) leading to a 
worse match between theory and experiments. 

The same methodology was used to determine the value of m in the non-specific correction 
for the free energy of formation of one base pair, where the elastic parameters are known but 
the free energy is unknown (see section 3.5 and Fig. 4 in the main paper): we look for the best 
matching between theory and experiments. 
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Figure 4: Determination of P. Data obtained at 1 mM [Mg 2+ ]. Gray lines are evaluations of 
the kinetic barrier using Kramers theory, eq. ( 30 1 , for different values of the persistence length 
(P = 0.5, 0.6, 0.7, . . . , 1.5 nm). Experimental data is analyzed using also different values of the 
persistence length. In this case, at P — 0.8 ± 0.1 nm we obtain the best overlapping between 
theory and experiments. 



41 



13 UV Absorbance experiments 



In order to determine the effect of salt on the stability of the hairpin we obtained the melting 
profile of the molecule using UV absorbance at 260 nm. The melting temperature was measured 
at 70±1°C using a buffer containing 100 mM Tris.HCl, 1 mM EDTA and no NaCl neither MgCl 2 . 
Results can be observed in Fig. [5]A We calculated the first derivative of the absorbance as a 
function of temperature (Fig. [5j3) and observed several maximums along the resulting profile (see 
arrows), which denote the presence of pre-melted states. For instance, regions with a richer A-U 
content in the middle of the stem may dissociate before the whole hairpin is unfolded (Fig. (5p). 
This result invalidates the two-states assumption used to extract thermodynamic parameters 
from the melting curve. 




T (°C) 

Figure 5: UV Absorbance of the RNA hairpin. (A) Experimental results of the absorbance 
of our RNA hairpin as a function of temperature. (B) First derivative of the absorbance as a 
function of temperature. Its maximum (black dashed line) defines the melting temperature. (C) 
The RNA hairpin under study has a region with a high A-U content in the stem which could 
lead to premelted states. 

We tried to obtain the melting profile adding 100 mM NaCl to the buffer, and the melting 
temperature was too high and the sample started boiling and evaporating before any relevant 
signal could be obtained. Therefore, for this RNA hairpin melting curves cannot be measured 
at the experimental conditions used in our pulling experiments with optical tweezers, and our 
results can only be compared with Mfold [15, 19, 68] and other theoretical predictions 81, 
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14 Rupture force histograms 
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Figure 6: Probability distributions of the unfolding and folding first rupture forces measured at 
different pulling speeds and different monovalent ionic condition. Red points are folding forces 
at 1.8 pN/s, green are unfolding forces at 1.8 pN/s, blue are folding forces at 12.5 pN/s and 
magenta are unfolding forces at 12.5 pN/s. 
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/ (pN) / ( P N) 

Figure 7: Probability distributions of the first rupture forces measured at different pulling speeds 
and different magnesium concentration. Red points are folding forces at 1.8 pN/s, green are 
unfolding forces at 1.8 pN/s, blue are folding forces at 12.5 pN/s and magenta are unfolding 
forces at 12.5 pN/s. 
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15 Tightly Bound Ion Model 

In this section the empirical equations of the Tightly Bound Ion (TBI, [SSI [55]) model used to 
predict the hairpin free energies at different ionic conditions are summarized [lU [45] . 
The external parameters, obtained from the Mfold server [T51 [THl EH] , are: 

Bases in the loop Ni=A 

Bases in the helix iV=20 

Hairpin diameter d=1.7 nm 

Interphosphate distance a=0.6 nm 

Enthalpy at 1 M [Mon+], M [Mg 2 +] A# =199 kcal/mol 

Entropy at 1 M [Mon+], M [Mg 2+ ] AS , =527.16 mkcal/Kmol 



In what follows x is the concentration of monovalent salt and y is the concentration of 
magnesium ions. Both parameters are given in units of M. Temperature T is given in Celsius. 
The empirical set of equations are: 



In Z l Mon (x) = a[(x) log(JV z - a/d + 1) + ^(^(JVj - a/d + l) 2 - b[{x) (33) 

lnZ c Mon (x) = c l 1 (x)N l -d[(x) (34) 

a\ (x) = (0.02Ni - 0.026) log (a;) + 0.54iV ; + 0.78 (35) 

b[( X ) = (- ^ + 0.006) log(x) - I ^ - 0.01 (36) 

c[(x) =0.07 log(ac) + 1.8 (37) 

4(a;)=0.211og(a:) + 1.5 (38) 

G l Mon ( x ) = - (l n Z Mon ( x ) ~ ln Z° Mon (x) ) (39) 

G h Mon (x,T) = H Q -(T + 273.15)5 Mo n(a:)0.001 (40) 

S M on(x) =S - 3.22(N - l) 9l (x) (41) 

g 1 (x) = a h 1 (x)+b h 1 (x)/N (42) 

a \{x) = -0.075 log(ac) + 0.012 log 2 (x) (43) 

b\{x) = 0.018 log 2 (a;) (44) 

Gmou(x, T) = G Mon (x, T) + G Mon (x) (45) 
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\nZ l Mg (y) = 4(y) log(7V z - a/d + 1) + 6£(y)(iV, - a/d + l) 2 - b l 2 (y) (46) 

\nZ c Mg {y)=c l 2 {y)N l -d l 2 {y) (47) 

4(y) = f- lvF 1 — + 0.32^) log(y) + 0.77V, + 0.43 (48) 



b l 2 {y) = 0.0002(JV ; + 1) log(y) - 5.9/{N t + l) 2 - 0.003 (49) 

4(y) = 0.0671og(y)+2.2 (50) 

4(t/) = 0.163 log(y) + 2.53 (51) 

G l Mg (y) = -(In Z l Mg (y) - lnZ c Mg (y)) (52) 

(53) 

G h Mg {y,T) =H - (T + 273.15)5 M9 (y)0.001 (54) 

S Mg (y) = S -3.22(N-l)g 2 (y) (55) 

g 2 (y) = a h 2 (y) + b h 2 (y)/N 2 (56) 

4(2/) = -0.6/N + 0.025 log(y) + 0.0068 log 2 (y) (57) 

=log(2/) + 0-38 log 2 (y) (58) 

GMg(y, T) = G h Mg (y, T) + G l Mg (y) (59) 

x' 1 (a;,i/)=!c/(a; + (7.2-20/JVi)(40-Iog(a:))») (60) 

C,m 9 (^!/) =^i(^y)G' Mo nW + ( 1 -4(^2/))G^(2/) (61) 

Gm„„, m >, y, T)=H -(T+ 273.15)S M on,Mg(x, y)0.001 (62) 

S Mon ,Mgfa y) = S - 3.22 [(iV - 1) (x^x, y) 9l (x) + (1 - a# (a, y))<&(y)) + <?i, 2 (z, 2/)] (63) 

Xkl{x ' V) = x + (8.1-32.4/iV)(5.2-log(x))y (64) 

3i, 2 (x, y) = -0.6a;f (x, y)(l - xfa, y)) \og{x) log((l/4 (a, y) - l)a;)/iV (65) 

GMon,Mg(x, y, T) = G , l IonMg (x, y, T) + G l Mon Mg (x, y) 



(66) 



Where GMon,M g {x,y 7 T) is the free energy at any temperature and at any monovalent and 
magnesium ion concentration. 
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16 Comparison to the counterion condensation theory 



There are two successful theories to account for the energetic interactions between ions in solution 
and nucleic acids: the Poisson-Boltzmann theory and the counterion condensation theory derived 
by Manning [HTl |H7j . These theories are based on different mean field approaches and neglect 
any kind of correlations between the ions in the solution. More recently a new theory known as 
the Tightly Bound Ion (TBI) model has been introduced 96 , which accounts for the different 
modes of correlations between counterfoils. 

In Fig. [8] we see the prediction provided by the Manning theory and the TBI model to the 
free energy of formation of our RNA hairpin as a function of the salt concentration. Because 
correlations between monovalent ions are negligible, we see that both the Manning theory and 
the TBI model give similar results under this condition (Fig. [8]A_). However, correlations between 
Mg 2+ are important and the TBI model gives an improved prediction in this case (Fig. [8)3) [56] . 



A B 




Figure 8: Comparison between experimental data, the counterion condensation theory and the 
TBI model to predict the behavior of the free energy of formation of our RNA hairpin as a 
function of salt concentration. 
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